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The viscous Cardassian cosmology is discussed, assuming that there is a bulk viscosity 
in the cosmic fluid. The dynamical analysis indicates that there exists a singular curve in 
the phase diagram of viscous Cardassian model. In the viscous PL model, the equation- 
of-state parameter is no longer a constant and it can cross the cosmological constant 
divide u>a = — 1, in contrast with same problem of the ordinary PL model. Other models 
possess with similar characteristics. For MP and exp models, evolves more near —1 
than the case without viscosity. The bulk viscosity also effect the virialization process of 
a collapse system in the universe: is increasingly large when the bulk viscosity is 

increasing. In other words, the bulk viscosity retards the progress of collapse system. In 
addition, we fit the viscous Cardassian models to current type la supernovae data and 
give the best fit value of the model parameters including the bulk viscosity coefficient r. 

Keywords: Bulk viscosity; Cardassian universe;cosmological dynamics; virialization. 



1. Introduction 

The current accelerating expansion of the universe indicated by the astronomical 
measurements from high redshift supernovae^ as well as accordance with other 
observations such as the Boomerang/Maxima/WMAP data^l and galaxy power 
spectra becomes one of the biggest puzzles in the research field of cosmology. One 
popular theoretical explanation approach is to assume that there exists a mysterious 
energy component, dubbed dark energy, with negative pressure, or equation of state 
with w = p/p < that currently dominates the dynamics of the universe(see^ and 
references therein). Such a component makes up 70% of the energy density of the 
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universe yet remains elusive in the context of general relativity and the standard 
model of particle physics. The most natural dark energy candidate is a cosmolog- 
ical constant which arises as the result of a combination of general relativity(GR) 
and quantum field theory. However, its theoretical value is between 60 — 120 orders 
of magnitude higher than the observed value. An alternative candidate is a slowly 
evolvingand spatially homogeneous scalar field, referred to as " quintessence" with 
w > — ]P and "phantom" with w < — ]P2 respectively. Since current observational 
constraint on the equation of state of dark energy lies in a relatively large around 
the so-called cosmological constant divide w\ — —1, it is still too early to rule out 
any of the above candidates. However, the expectation of explaining cosmological 
observations without requiring new energy component is undoubtedly worth of in- 
vestigation. GR is very well examined in the solar system, in observation of the 
period of the binary pulsar, and in the early universe, via primordial nucleosynthe- 
sis, but no one has so far tested in the ultra-large length scales and low curvatures 
characteristic of the Hubble radius today. Therefore, it is a priori believable that 
Friedmann equation is modified in the very far infrared, in such a way that the 
universe begins to accelerate at late time. Freese and Lewis® construct so-called 
Cardassian universe models that incarnates this hope. The Cardassian universe is a 
proposed modification to the Friedmann equation in which the universe is flat and 
accelerating, and yet contains only matter(baryonic or not) and radiation. But the 
ordinary Friedmann equation governing the expansion of the universe is modified 
to be 



where p consists only of matter and radiation, H is the Hubble "parameter" which 
is a function of time, a is the scale factor of the universe, and G = l/m^ = ^ is 
the Newtonian gravitational constant. Note that as required by inflation scenario 
and observations of CMB, the geometry of the universe is flat, therefore, there are 
no curvature terms in the above equation. The function g{p) returns to ordinary 
Friedmann equation at early times, but that takes a distinct form that drives an 
accelerated expansion in the recent past of the universe at z < 0(1). In Cardassian 
models, we simply set the cosmological constant A = 0, and the only compositions 
are matter and radiation. A possible interpretation of Cardassian arose from consid- 
eration of braneworld scenarios, in which our observable universe is a 3-dimensional 
brane embedded in a higher dimensional universe. An alternative interpretation 
had been also discussed^, in which one developed a 4-dimensional fluid descrip- 
tion of Cardassian cosmology. The observational constraints of Cardassian models 
have been extensively studied^. The simplest model is power law Cardassian mo- 
dle(PL) with an additional term p n , which satisfies many observational constrains: 
the first Doppler peak of CMB is slightly shifted, the universe is rather older, and 
the early structure formation(z > 1) is unaffected. Furthermore, one can consider 
other forms of g{p) including the modified polytropic model(MP) and the exponen- 
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where p car d is a characteristic constant energy density and q and n are two di- 
mensionless positive constants^!. Obviously, at early times, p is much larger than 
the characteristic energy density p car d, g{p) — * P, i.e. Eq.(JTJ) recovers the ordinary 
Fricdmann equation. 

On the other hand, the dissipative effects, including both bulk and shear vis- 
cosity, are supposed to play a very important role in the astrophysics (see Ref llll 
and references therein) and the nuclear physics (see Ref 1T21 and references therein). 
Under the conditions of spatial homogeneity and isotropy, a scalar bulk viscous 
pressure is the solely admissible dissipative phenomenon. The viscosity theory of 
relativistic fluids was first suggested by EckarlP^l and Landau and Lifshit J^, wn0 
considered only first-order deviation from equilibrium, which leads to parabolic dif- 
ferential equations and hence to an infinite speed of propagation for heat flow and 
viscosity, in contradiction with the principle of causality. The relativistic second- 
order theory was founded by Israel^, and has also been used in the evolution of the 
univers However, the character of the evolution equation is very complicated 
in the framework of the full causal theory. Therefore, the conventional theorjEl is 
still applied to phenomena which are quasi-stationary, i.e., slowly varying on space 
and time scales characterized by the mean free path and the mean collision time. In 
the case of isotropic and homogeneous cosmologies, the dissipative process can be 
modelled as a bulk viscosity Q within a thermodynamical approach. Some original 
works on the bulk viscous cosmology were done by Belinsky and Khalatnikov^l 
The bulk viscosity introduces dissipation by only redefining the effective pressure, 
Peff, according to p e ff = P — 3£H where H is the Hubble parameter. The condi- 
tion £ > assures a positive entropy production in conformity with the second law 
of thermodynamics^. We are interested in two solvable cases: (i) £ = \ /r 3n~ 1 TH , 
where r is a constant. This assumption implicates that £ is directly proportional 
to the divergence of the cosmic fluid's velocity vector. Therefore, it is physically 
natural, and has been considered previously in an astrophysical contexlM^ (ii) 
C = T (g(p)) a+ ^ ■ This dependence is more complicated, but one can see in the 
following some interesting results obtained. Obviously, case (i) is equivalent to one 
with a = — 1 of case (ii). 

In this paper, we consider a viscous Cardassian cosmological model for the ex- 
panding universe, assuming that there is bulk viscosity in the cosmic fluid. The 
dynamical analysis indicates that there exists a singular curve in the phase diagram 
for various models. The numerical result show that is increasingly large when 
the bulk viscosity is increasing. Furthermore, we fit the viscous Cardassian models 



19, 2009 0:5 WSPC/INSTRUCTION FILE vcu 



4 CHANG-BO SUN,JIA-LING WANG and XIN-ZHOU LI 

to current SNela data and give the best fit values of the parameters including the 
bulk viscosity coefficient r. 

2. The model 

In the viscous Cardassian model, the modified Friedmann equation is described by 
Eq.(p}. For the late-time evolution of the universe we neglect the contribution of 
radiation, so that Eq.([T]) is reduced to 

H 2 = y 5 (p m ) (3) 

where p m is the energy density of matter, which keeps conserved during the expan- 
sion of the universe, i.e., 

p m + 3H(p m - 3#Cm) = (4) 

where £ m is the bulk viscosity for the matter energy density p m . In the £ m = case, 
the evolution of matter takes the ordinary manner p m = p mt o(l + z) 3 where p m ,o 
is the present value of energy density of matter. However, the evolution of matter 
is not likely power law of (1 + z) in the £ m ^ case. Similarly, the conservation 
equation of total energy density can be written as 

g(p m )+3H[g(p m )+p-3H(}=0 (5) 

where £ is the bulk viscosity for total energy density g{p m ). Combing Eqs.(|l|) and 

dpr, 



<j5j> , it is easy to check that £ = d9 ^ pm ) C, m ■ Following RefjH we take the energy 



density to be the sum of two terms: 

g(p m )=Pm+Pk (6) 

where pk is a Cardassian contribution. The thermodynamics of an adiabatically 
expanding universe tell us that pressure of Cardassian contribution is 

dg(pm) , v 

Vh = Pm—^ g(Pm), (7) 

and p = pk because of p m = 0. From Eqs.© and 0, we have the equation of state 

Pk - 3ffC rjn 
w k = —, — ^ ■ (8) 

9{Pm)- Pm 

Furthermore, we have 

2 

H + H 2 = - = — [g( Pm )+Z{Pk-ZHQ] (9) 
a 6 

which is the acceleration equation of the cosmological expansion. Using Eqs.([3|),([5|) 
and ((?]), we obtain 

dg(p m ) . „ dg(p m ) r— 

a — h3p m — 3KTj3g(p m K = °- ( 10 ) 

da ap rn 
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We shall be interested in the evolution of the late universe, from t = to onwards, 
where to is the initial time and the corresponding scale factor a(to) — ao and energy 
density p m (to) = Pm,o- From Eq. lfTO]) . we have 



a = ao exp 



dg{p m ) 



Pm.o 3K^3g(p m )( - 3p m d9 }p™ ] 



(11) 



which is the general relation between the cosmological scale factor a and the energy 
density p m . 

In what follows, we focus on the viscous Cardassian model, as the cosmological 
dynamics are analytically solvable when we choose ( — \/3k~ 1 tH where r is a 
constant. The scale factor is given by 

a _ /(O m ) 



ro.OJ 



(12) 



where Q m is the density parameter of matter and £l m ,o is the present value of il ri 
For the PL model, we have 



/(n. 



(V3KT-1 



3(»-l)(r,-^r)(/3KT-l) 



^ [n{n m - 1) - O m + V3KT]^(n-l 

and the equation of state 

(1 — n)(Q m — 1) — \fSnT 



i - a 



(13) 



(14) 



In the late time of universe, the new term of PL model is so large that the ordinary 
first term p m can be neglected. In other words, we have f2 TO < 1, so that the 
expansion is superluminal (accelerated) for n < | + \/3kt. 
For the MP model, we have 



f(M = 



9. 



(9.? n - 1)»+\/3kt-1 



and 



(V3kt + n - 1 - ntt q m )^ KTn 

n{W m - 1) - V3kt 



Wk 



For the exp model, we have 

/(Am) 



1-0 T , 



ln£l m 



and 



(1- V3kt + nlnSl m )^ K 
n In fl m — \^3kt 



3n(l- s/3kt) 



w k = 



i - a 



(15) 



(16) 



(17) 



(18) 



We plot the evolution of Wk for different values of kt in Fig. [T] As can be seen, 
in viscous PL model, Wk is no longer a constant and it is dependent on time that 
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can cross the cosmological constant divide w\ = — 1. For MP and exp models, Wk 
evolves more near —1 than the case without viscosity. 

In the r = case, Eqs. (fT3l) - (IT8t reduce to the results of ordinary Cardassian 
model without viscosity and Eq. (fT2| is solvable for the three models: 



;(!+*) 



3(l-n) 



l+T^(l + -) 3<1 - ) ' 

1 "m,0 

l+-^-(l + ^) 3<!n 



for PL 



for MP 



(19) 



(1+*)- 

m,0 



for exp 



Fig. 1. In the £ = VSk^tH case, the evolution of the equation-of-state parameters uij. for 
different values of nr. 



3. Autonomous system 

A general study of the phase space system of quintessence and phantom in FRW uni- 
verse has been given in Refl20l For the viscous Cardassian cosmological dynamical 
system, the corresponding equations of motion can be written as 

-»"<)• « 

Pm + 3H (p m - 3HC-^-) = (21) 



'dg(p m ), 

and Eq.Q. To analyze the dynamical system, we rewrite the equations with the 

« 2 Pm „ _ «fC _ _ dg(p, 
3H 2 >» ~ 3H ~ dp„ 



following dimensionless variables: x = \ ,y = ^rS,z — d9 j Pm ^ and N = In a. The 



dynamical system can be reduced to 

dx 9y „ 9 „ 

— = — -3x + 3x 2 z-9xy, 22 
aN z 

d V , 3 9 2 , OQ N 

"T77 = U + — — y , (2d) 

dTV 2 y 2 y ' v ' 

where z and u are both functions of x and y. z(x,y) is determined by the specified 
form of g(p m ) and u = jffjr is determined by the specified form of viscosity The 
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equation of state can be expressed in terms of the new variables as 

7 



xz-Zy-l 

W k = : (24) 



and the sound speed is 



for PL 



n(n — l)(l — x) 

c 2 = J »(l-^)[n-l + ( 9 -l)n^] f Mp (25) 



nx q + l — n 
n(l— n) In x-\-n 2 (In x) A 
1+n In a; 



for exp. 



We choose C = T (.9(/°)) Q 5 an d show the critical points of the autonomous 
systems and their properties in Table [TJ In Fig. [21 the phase diagrams of the three 
models are given, respectively. It is worth noting that there exists a singular curve 
(bold line) in the viscous Cardassian's phase diagrams. 



Table 1. The critical points and their properties in the £ = i~(g(pm)) 2 case. 



Critical points(a;,y) 


Eigenvalues 

PL model 








Stability 


(0,0) 
(1,0) 

(*o.£=iP»±S) ( 


(3(n-l),-3n(l + a)) 
(3(l-u),-3(l + a)) 

^0,3[(l- n )(2 +Q )x + (l+o)„-l +ri+(1 ^ ti)xo 


) i 


f_- 



or? 


{ 

^TH ^ 2 ^ 

l rT <n<^,Q 


0<n<| ,a>-l 

J , stable 
or n<0,a< — 1 

unstable 

, [l-4T 1 + 3n 2 -2(l-n) 2 x ]x -n 2 


MP model 


(1,0) 


(-3(l + (*),3<jn) 
(0,3q(n-l)-3) 

^0,3^(l + c)(l- re ) + (l + <j + Q ) n x«- T -^a_^J 


) 1 


(l+c)(l- 
and l - 


unstable 
J q(n — 1)<1, stable 
1 q(n— 1)>1, unstable 

i + iix«^0 


exp model 


(1,0) 
(e"n,0) 

(*o, 1+ " 3 ta *°) 


(3n,-3(l + a)) 
(0,3(n-l)) 

(o,s[l + a + n+n(l + a) lnx - 1 + „" n3!0 ] j 


unstable 
J 0<n<l, stable 
1^ n>l, unstable 

r 1+a+n+n(1+a)ln:co < TT _n_ ^ 

and 1 + 



4. Virialization 

The spherical collapse formalism developed by Gunn and Gott^ is a simple but 
powerful tool for studying the growth of inhomogeneities and bound systems in the 
universe. It describes how an initial inhomogeneity decouples from the expansion of 
the universe and then expands slower, eventually reaches the state of turnaround and 
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Fig. 2. In the f = t (g(p)) a+ 2 case, the phase diagrams of (x,y) for PL, MP and exp models. 
There is a singular curve (bold line) which formed of critical points in the phase diagrams. We 
choose the parameters to be n = —0.01, a = 0.2 for PL, q = 1.15, n = 1.1, a = 1.5 for MP and 
n = 0.8, a = 2.5 for exp, respectively. 



collapses. Physically we assumes that the collapse system goes through a virializa- 
tion process and stabilizes at a finite size. An important parameter of the spherical 

collapse model is the ratio between the virialized radius and the turnaround radius 

n mr 22 



Rta 



Using energy conservation between virialization and turnaround, one can get 



[U + T] mr = U ta 



(26) 



where U is the potential energy, T = ^ is the kinetic energy of the system. If 
the collapse object is made up of only one energy component, the potential and the 
potential energy within it are 



$(r) 



-2irG(l + 3w)p [R 



(27) 



and 



U = - I p^dV 



(28) 



where R is the radius of the spherical collapse system. In the Einstein-de Sitter 
universe, for a spherical perturbation with conserved mass M, <j>(r) = ^|(r 2 — 3i? 2 ), 



U = — gjjT , T V i r = and the ratio of virialization to turnaround is -^F 
In the PL model without bulk viscosity, we have 



$ = ~2GTi Pm \R 2 -^r 



and 



Pcard 



1 + 3 



n — 1 



(1 - n)p„ 



J m I year 



(29) 
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Pm ^~ Pcard 



Pcard 



Prn + (3n - 2)p card 



Pm 
Pcard 



(30) 
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Using the relation p m . V ir = Pm.taS, 3 where £ = Eq. (|26|) can be written as 

{2£ [1 + (1 + z ta ) 3 "- 3 (tt m | - 1)] [(3n - 2)(1 + z ta ) 3 - 3 (^ - 1) + 1] - 1} e n + 

(9n 2 - 15n + 4)(1 + z ta ) 3 ^ 3 (fi^ - 1)£ 3 + 3 " + (3n - 2)(6n - 7)(1 + z ta ) 6 ^ 6 (fi m | - l) 2 f 

= 0. (31) 

It is difficult to find analytical solution for Eq. (f3"Tj) , therefore we have to investigate 
it numerically. The result is shown in Fig. [3] (solid line). In the viscous PL model 
with £ = v / 3K _1 rif, same problem becomes more complex. As the function of scale 
factor a, p m is determined by 

3(v / 3«T-l)(v / 3nT-re) 




Pcard \ I -, / Pc 



Pm,0 



—n 


( 
















Pcard 


r 


Pm,0 





V3k 



K'3,vr-„)(^) (^J ^.(32) 

Combing (|32)) and (|26|) . we obtain the relation between £ and turnaround redshift 
Zta- In Fig- El we show f as a function of Zt a for different bulk viscosity coefficients. 
As can be seen, in Cardassian models, the ratio is always larger than i, and 
it get larger and larger with the evolution of the universe. It means that collapse 
process will be harder and harder to occur (£ — > 1). Furthermore, £ becomes larger 
when r takes larger value. 



Fig. 3. In the viscous PL model with £ = V3k~ 1 tH, £ as a function of zta for different bulk 
viscosity coefficients r. 



5. Fit the model parameters to supernovae data 

In general, the approach towards determining the expansion history H(z) is to as- 
sume an arbitrary ansatz for H{z) which is not necessarily physically motivated 
but is specially designed to give a good fit to the data for cIl(z). Given a particular 
cosmological model for H(z; ax, a n ) where a\,...,a n are model parameters, the 
maximum likelihood technique can be used to determine the best fit values of pa- 
rameters as well as the goodness of the fit of the model to the data. The technique 
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can be summarized as follows: The observational data consist of N apparent mag- 
nitudes mi(zi) and redshifts z, with their corresponding errors a mi and a Zi . These 
errors are assumed to be gaussian and uncorrelated. Each apparent magnitude mi 
is related to the corresponding luminosity distance di by 



m{z) = M + 5 \og 1 



d L (z) 



Mpc 



25, (33) 



where M is the absolute magnitude. For the distant SNela, one can directly observe 
their apparent magnitude m and redshift z, because the absolute magnitude M of 
them is assumed to be constant, i.e., the supernovae are standard candles. Obviously, 
the luminosity distance di,{z) is the 'meeting point' between the observed apparent 
magnitude m(z) and the theoretical prediction H(z). Usually, one define distance 
modulus (i(z) = m(z) — M and express it in terms of the dimensionless 'Hubble- 
constant free' luminosity distance Dl defined hyD^z) = HocIl(z) / c as 

M(z) =51og 10 (£>i(2)) + /x , (34) 

where the zero offset /io depends on Hq (or h) as 

Mo = 5 log 10 (^7) + 25 = -5 log 10 h + 42.38. (35) 

The theoretically predicted value D t ^(z) in the context of a given model 
H(z; a%, a n ) can be described by 

D?(z) = {l + z) f dz' -— ^ -. (36) 

Jo H{z';a 1 ,...a n ) 

If we assume prior to the parameters Q m ,o , the viscous Cardassian models predict 
a specific form of the Hubble parameter H (z) as a function of redshift z in terms 
of two parameters n and r for the viscous PL model and exp model. Therefore, the 
best fit values for the parameters (n, r) are found by minimizing the quantity 

2 , v y> K*(^) - 51og 10 £>ff(z«; n, r) - 2 
X (n, t) = -3 ■ (37) 

i=l » 

Since the nuisance parameter fi is model-independent, its value from a specific good 



1001 

fit can be used as consistency test of the data and one can choose a priori value 
of it (equivalently, the value of dimensionless Hubble parameter h) or marginalize 
over it thus obtaining 

X 2 (n,r) = A(n,r) - +ln (£\ , (38) 



where 



A (n, t) = -2 — , (39) 
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N 



[^(z. t )-51og 10 ^ 



)] 



(40) 



and 




(41) 



In the latter approach, instead of minimizing x 2 (n,r), one can minimize x 2 (n,r) 
which is independent of /irj. 

We now apply the above described maximum likelihood method for the viscous 
PL model using Gold dataset which is one of the reliable publish dataset consisting 
of 206 SNela (N=206j231. In Fig.Hl contours with 68.3%,95.4% and 99.7% confidence 
level are plotted, in which we take a marginalization over the model-independent 
parameter fiQ. The best fit as showed in the figure corresponds to n = —0.252 and 
kt = 0.02, and the minimum value of x 2 is 169.739. 

In Fig.0 contours with 68. 3%, 95. 4% and 99.7% confidence level are plotted for 
the viscous exp model. The best fit is n — 0.8 and kt = 0.074, and Xmin ~ 169.857. 
Obviously, the allowed ranged of the parameters n and r favor that there exists an 
effective phantom energy in the universe. 

As for the viscous MP model, using the above marginalization method, we find 
the minimum value of % 2 (g, n,r) is 169.823 and the corresponding best fit value of 
parameters are q — 0.233, n — 3.853, and kt = 0.06. 



Fig. 4. The 68. 3%, 95. 4% and 99.7% confidence contours of parameters n and kt using the Gold 
SNela dataset and marginalizing over the model-independent parameter no with a prior £l m ,o = 0.3 
in PL model. 



6. Discussion and conclusion 

In above sections, we have investigated the viscous Cardassian cosmology for the 
expanding universe, assuming that there is a bulk viscosity in the cosmic fluid. We 
consider two solvable cases: (i) ( = v / 3k~ 1 tH; (ii) £ = t {g{p)) a+ " s . In case (i), 
there exist exact solutions of a as a function of fi m ,o for PL, MP and exp models. 
Contrary to the naive PL model, in viscous PL model, the effective equation-of-state 




1 



Pi 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



Fig. 5. The 68. 3%, 95. 4% and 99.7% confidence contours of parameters n and kt using the Gold 
SNela dataset and marginalizing over the model-independent parameter po with a prior Q m ,o = 0.3 
in cxp model, w^ q and w t fi are equation-of-state parameters of p^ and g(p m ) at present time, 
respectively. 

parameter Wk is no longer a constant and it is dependent on time that can cross the 
cosmological constant divide wa = — 1 from Wk > — 1 to Wk < — 1. Other models 
possess with similar characteristics. For MP and exp models, w fc evolves more near 
— 1 than the case without viscosity. Moreover, the dynamical analysis indicates that 
there exists a singular curve in the phase diagram of viscous Cardassian . 

On the other hand, the bulk viscosity can effect the collapse process of a bound 
system in the universe. We consider the virialization of a spherical collapse model 
using the spherical collapse formalism in PL Cardassian. The numerical result in- 
dicates that the parameter £ = is always larger than \ and get larger and 
larger with the evolution of the universe. Furthermore, £ becomes larger when r 
is increasing that means the bulk viscosity retards the progress of collapse system. 
Obviously, the bulk viscosity should be small to insure that the viscous cosmology 
theory isn't in contradiction with the galaxy formation theory. 

Using maximum likelihood technique, we constrain the parameters of viscous 
Cardassian models from the supernova data. If we assume prior that £l m .o = 0.3 
as indicated by the observation about mass function of galaxies, the best fits of kt 
are 0.02, 0.06 and 0.07 for PL, MP and exp models, respectively. In the following 
works, we plan to use other cosmological and astrophysical observations such as 
CMB, BAO and LSS to further constrain the viscous Cardassian parameters and 
the bulk viscosity coefficient. 
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